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A systematic investigation of the nuclear-polarization effects in one- and few-electron heavy ions is presented. 
The nuclear-polarization corrections in the zeroth and first orders in 1/Z are evaluated to the binding energies, 
the hyperfine splitting, and the bound-electron g factor. It is shown, that the nuclear-polarization contributions 
can be substantially canceled simultaneously with the rigid nuclear corrections. This allows for new prospects 
for probing the QED effects in strong electromagnetic field and the determination of fundamental constants. 
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The enormous progress made in experimental investiga¬ 
tions of heavy high ly charged ions during the last decades 
(see, e.g., Refs. Uj-lst] and references therein) has triggered the 
vigorous development of ab initio QED theory in the presence 
of strong nuclear fields. The relativistic behaviour of electrons 
in highly charged ions requires a fully relativistic description 
from the very beginning, i.e. nonperturbative in the aZ pa¬ 
rameter, where Z is the nuclear charge number. This plays a 
key role in contrast to QED theory for light atomic systems, 
where the parameter aZ is employed as an expansion param¬ 
eter. Over the last decades an essential progress has been 
achieved in theoretical calculations of various spectroscopic 
properties of highly charged ions, such as transition energies, 
hyperfine splitting (HES), and g factor (see Refs. |@-@1 for re¬ 
views). In many cases further improvement of the achieved 
theoretical accuracy seems strongly limited by the lack of the 
knowledge of the nuclear properties. E.g., in the case of the 
g factor of H-like lead ion the uncertainty of nuclear charge 
distribution correction is the main source of the total uncer¬ 
tainty, and in the case of the HES in the H-like bismuth ion 
the uncertainty of the nuclear magnetization distribution cor¬ 
rection (so-called Bohr-Weisskopf effect) strongly masks the 
QED contributions. In Ref. it was proposed to consider a 
specific difference of the HES values of H- and Li-like ions 
with the same nucleus, where the uncertainty of the Bohr- 
Weisskopf effect is significantly reduced and the QED effects 
can be tested on the level of a few percent. In the case of the 
g factor similar cancellations of the finite nuclear size correc¬ 
tions have been recognized for the specific differences of the 
g factors of H- and Li-like ions in Ref. ifl^ and of H- and B- 
like ions in Ref. fil , respectively. These differences can be 
calculated with a substantially higher accuracy, which opens 
excellent perspectives for a test of the QED effects and even 
provides a possibility for an independent determination of the 
fine structure constant from the strong-held QED theory. 

Another nuclear effect appears due to the intrinsic nu¬ 
clear dynamics, where the nucleus interacting with electrons 
via the radiation held can undergo real or virtual electro¬ 
magnetic excitations. The latter effect leads to the nuclear- 



FIG. 1: Feynman diagrams representing the lowest-order nuclear- 
polarization effect to the electron binding energy. The bound electron 
(double line) interacts with the nucleus (heavy line) in its ground 
state via the exchange of virtual photons (wavy lines). 


polarization (NP) correction, e.g., to the binding energy of the 
electrons. Being restricted to phenomenological descriptions 
of the nucleon-nucleon interaction the NP correction sets the 
ultimate accuracy limit up to which QED corrections can be 
tested in highly charged ions. Therefore, an important ques¬ 
tion should be addressed: To which extent can NP corrections 
be canceled in specific differences? In this Letter, we rigor¬ 
ously examined the NP and the screened NP cotTections to the 
binding energies, HES, and g factor of heavy highly charged 
ions. We analyze the ratio of the finite nuclear size and the NP 
corrections and consequently evaluate the NP contribution to 
the specific differences, designed for the cancellation of the 
finite nuclear size and the Bohr-Weisskopf effect. 


In Refs. ||12lll3|] a relativistic field theoretical approach to 
the NP corrections in electronic atoms incorporating the ef¬ 
fects due to virtual collective nuclear excitations within the 
framework of bound-state QED for atomic electrons was de¬ 
veloped. This approach was successfully applied in calcula¬ 
tions of NP corrections to the binding energies EMIgo the 
HES d, and to the bound-electron g factor llisll . 


The lowest order diagrams describing the NP effect are de¬ 
picted in Pig.[Tl As virtual nuclear excitations we account for 
the dominant ones arising from the collective nuclear dynam¬ 
ics, such as rotations of deformed nuclei, harmonic surface 
vibrations, and giant resonances. Since the velocities asso¬ 
ciated with the collective nuclear dynamics are nonrelativis- 
tic we can restrict to the nuclear charge-density fluctuation 
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(electric multipole transitions) and neglect contributions aris¬ 
ing from fluctuations of the nuclear vector current (magnetic 
multipole transitions). It is most suitable to employ Coulomb 
gauge and keeping the longitudinal component Dqq of effec¬ 
tive photon propagator only. It describes the interaction 
between the electrons and virtual nuclear transitions and takes 


the form Ill2l 11311 : 


Z)oo(i’i,i'2,a;) = y^B{EL;L 
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Here wl is the nuclear excitation energy and B{FL\ L —>• 0) 
is corresponding reduced electric transition probability. This 
form of the propagator is also very suitable for numerical cal¬ 
culations, since it exclusively depends on phenomenological 
quantities such as transition energies ujl, and corresponding 
electric transition strengths. The radial dependence carried by 
the functions Fl may be specified utilizing, e.g., a sharp sur¬ 
face model for describing the corresponding collective nuclear 
multipole transition densities Ii3-i6|. The nuclear ground- 
state sphere radius Rq is determined as Rq — i/b/S 
where is the root-mean-square charge radius. 

The energy shift due to the lowest-order NP effect is given 
by 


Aii/Np = 


27r 
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where the summation runs over the complete Dirac spec¬ 
trum, and u = 1 — iO preserves the proper treatment of 
poles of the electron propagator. In Table |T] the leading or¬ 
der NP corrections are presented for the Is, 2s, and 2pi/2 
binding energies in | 2 ®Pb and g|®U ions. For the low-lying 
rotational and vibrational nuclear excitations the experimen¬ 
tal values for the excitation energies ujl and electric transi¬ 
tion strengths B{FL]L 0) ^e taken from Ref. iS for 
the ion and from Ref. for the ion. The cor¬ 

responding wl and B{FL\ L ^ 0) values for the giant res¬ 
onances have been estimated employing the phenomenologi¬ 
cal energy-weighted sum rules 1211] . The summation over the 
spectrum of the Dirac equation has been performed employ¬ 
ing the dual-kinetic-balance finite basis set method 12211 with 
basis functions constructed from B splines lEsIl . The Dirac 
spectrum is calculated in the field of extended nuclei utilizing 
nuclear charge density distributions with recent values for the 
radii = 5.5010 fm and = 5.8569 fm in the 

case of lead and uranium ions, respectively. As one can see 
from Table U the obtained results are in a fai r ag reement with 
the previous calculations presented in Ref. Jl^, and a better 
agreement is found with the results obtained by the direct nu¬ 
merical integration. 

In Table U we also present the corresponding leading order 
finite nuclear size corrections AiJpNS together with a ratio 
of the NP and finite nuclear size terms A^p /fns defined as 

^NP/FNS = 


NP 


/Ai^FNS- These ratios appear to behave 


TABLE I: Nuclear-polarization ALInp and finite nuclear size 
ASfns corrections to the Is, 2 s, and 2 pi /2 binding energies in 
g 2 Pb and 92 U ions. The nuclear-polarization corrections are com¬ 
pared with the previous calculations (H. The corresponding ra¬ 
tios of the nuclear-polarization and finite nuclear size contributions 
Anp /fns ^6 presented. 


State 

Is 

2s 

2Pl/2 


82 ® Pb nucleus 



Ai^NP (meV) 

-28.89 

-5.033 

-0.4249 


-29.3“ 

-5.0“ 



-STS® 

- 5 . 5 '’ 


ASfns (eV) 

67.18 

11.66 

0.9991 

Anp/fns (10 

■®) -0.430 

-0.431 

-0.425 


92 ®U nucleus 



ASnp (meV) 

-188.2 

-35.88 

-4.153 


-197.6“ 

-37.2“ 

-4.2“ 


- 213 . 4 '’ 

-40.9*’ 

-4.6*’ 

Ai?FNS (eV) 

198.6 

37.73 

4.412 

Anp/fns (10 

■3) -0.947 

-0.951 

-0.941 


“Ref. |T6J: direct numerical integration (a). 
'’Ref. d: B -spline calculations (b). 





FIG. 2: Feynman diagrams representing the screened nuclear- 
polarization effect to the electron level energy. Notations are the 
same as in Fig.[T] 


rather similar for all the considered electron state. This means 
that canceling the finite nuclear size corrections in an energy 
difference the NP effect will be also reduced to a large extent. 

However, the hydrogenic excited energy states are not al¬ 
ways accessible experimentally, e.g., at present, the highest 
accuracy was achieved in measurements of the 2 ni /2 ~ 2s 
transition energy in heavy Li-like 238pj89-i- j24|] . There¬ 

fore, it becomes of distinct importance to investigate the NP 
effect in few-electron ions. In the presence of other elec¬ 
trons, in addition to the leading order one-electron NP cor¬ 
rection, terms combining the interelectronic-interaction and 
NP effects appear. In analogy to the corresponding QED cor¬ 
rections we refer to them as the screened NP contribution. 
The diagrams of the NP corrections to the one-photon ex¬ 
change are depicted in Fig. |2] Expressions for the energy 
shift due to this effect can be derived according to the Eeyn- 
man diagrams depicted but are to lengthy to be presented 
here. Their evaluations have been performed in both Eeynman 
and Coulomb gauges for the photon propagator describing 
the electron-electron interaction, thus providing an accurate 
check of the numerical procedure. The results obtained for the 
screened NP corrections to the (Is)^ binding energy and to 
the (ls)^2s, (ls)^ 2 pi/ 2 , and (ls)^ 2 p 3/2 ionization energies 
in g 2 ®Pb and gf®U ions are presented in Table HJ As one can 
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TABLE II: Screened nuclear-polarization and screened finite nuclear 
size corrections A_Esnp and ASsfns, respectively, to the (Is)^ 
binding energy and to the (ls)^ 2 s, (ls)^ 2 pi/ 2 > snd (ls)^ 2 p 3/2 ion- 
ization energies (with the opposite sign) in 32 Pb and 92 U ions. The 
ratio of the screened nuclear-polarization and screened finite nuclear 
size contributions Asnp/sfns are presented. 


State 


(ls)^2s 

(ls)"2pi/2 

(ls)"2p3/2 


82® Pb nucleus 



ASsnp (meV) 

0.8441 

0.3017 

0.1218 

0.0335 

ASsfns (eV) 

-1.9629 

-0.6988 

-0.2763 

-0.0863 

Asnp/sfns (10 ®) 

-0.430 

-0.432 

-0.441 

-0.388 


238i 
92 ' 

U nucleus 



ASsnp (meV) 

5.498 

2.048 

0.9369 

0.1793 

ASsfns (eV) 

-5.802 

-2.152 

-0.9816 

-0.2117 

Asnp/sfns (10 ^) 

-0.948 

-0.952 

-0.955 

-0.847 


see, the screened NP correction is comparable with the leading 
order term, and has to be taken into account in specific differ¬ 
ences constructed for eliminating the finite nuclear size effect. 
In Table im we also present the finite nuclear size effect com¬ 
ing from the first-order interelectronic-interaction correction, 
the so-called screened finite nuclear size correction AT^gpNS^ 
together with the corresponding ratio Ag^p /sfns defined as 
Asnp/SFNS = ^^'Snp/A^'Sfns- The ratio of screened NP 
and finite nuclear size corrections appears to be rather stable 
for different electronic configurations. This opens a possibil¬ 
ity to eliminate in such differences not only the finite nuclear 
size corrections, but also the NP terms to a rather large extent. 

As an example for such a cancellation, let us consider the 
following difference. One of the most precise measurements 
were performed for the Is Lamb shift in H-like ura¬ 
nium and for the 2 pi /2 — 2s transition energy 

Aii;( 2 pi/ 2 - 2 s) jjj u89-i- 124], In both cases the uncertainty of 
the finite nuclear size correction essentially contributes to the 
total theoretical error bars. Thus, we can construct the follow¬ 
ing specific difference A'E — Ai?(^PL 2 - 2 s) _|_ ~ 

355.8 eV, where the parameter ^ = 0.161856 is chosen in a 
way to cancel the leading order and the screened finite nuclear 
size terms. Such a cancellation is also rather stable with re¬ 
spect to the employed nuclear charge distribution model. The 
NP corrections are canceled in A'E up to 10“^ eV opening 
a possibility for unprecedented tests of strong-held QED. Al¬ 
though we have considered here only nuclear excitations of 
electric type, this conclusion will hold in general, since the 
cancellation being discussed is observed for each individual 
nuclear excitation. In view of this, we can expect, that this 
cancellation will be rather independent from the employed nu¬ 
clear models. 

Let us now turn to the NP effects to the HFS in few-electron 
ions. The HFS transition line in Li-like Bi®°+ has been re¬ 
cently observed and directly measured in a laser spectroscopy 
measurement at GSI 10.0]. This has allowed for the first time 
to compare experimental and theoretical values for the spe¬ 
cific difference between HFS of H- and Li-like bismuth ions. 








FIG. 3: Feynman diagrams representing the lowest-order nuclear- 
polarization effect in the presence of an external potential. The 
dashed line terminated with the triangle denotes the interaction with 
the external magnetic field. Notations are the same as in Fig.[T] 



FIG. 4: Feynman diagrams representing the screened nuclear- 
polarization effect in the presence of an external potential. For 
brevity, we depict here only the direct part of the nuclear-polarization 
correction for all the diagrams except the first one. Notations are the 
same as in Fig.[T] 


Although the present experimental accuracy is smaller than 
the theoretical one 126, 27] in the future SPECTRAP Penning 
trap facility this will be improved by three orders of magnitude 
i28ll . A substantial progress in the theoretical calculations of 
the specific difference |2^ 27, ^ 33] allows to improve the 
theoretical accuracy by an order of magnitude, and the present 
uncertainty is partially restricted by the NP correction, which 
was so far known only for the Is HFS 

In this Letter we present for the first time results for the 
leading order and the screened NP corrections to the HFS 
of H-, Li-, and B-like bismuth ions. The leading order and 
the screened NP effect are given by the diagrams depicted in 
Figs.|3]and|4l respectively. 

The basic expressions for the leading order diagrams are 
similar to those derived in Ref. lfT3l . and for the screened 
diagrams they are rather bulky and will be presented else¬ 
where. The numerical procedure has been accurately checked 
by employing the Feynman and Coulomb gauges for the pho¬ 
ton propagator describing the interelectronic interaction. For 
the nuclear parameters of low-lying vibrational levels of the 
nearly spherical odd-even | 3 ®Bi nucleus we have employed 
the corresponding collective vibration levels in the neigh¬ 
boring even-even isotope of 82 ®Pb (weak coupling limit). 
Moreover, we have evaluated the effect of single-nucleon 
excitations and the effect going beyond the weak-coupling 
limit. Both have been found to be negligible in the case of 
Ig^Bi nucleus compared to the effect of collective core ex¬ 
citations. The detailed consideration will be presented else¬ 
where. The obtained results for the leading order and the 

















































4 


TABLE III: Individual contributions to the leading order AShfs.np 
and the screened AShfs.snp nuclear-polarization corrections to the 
ground state hyperfine splitting of H-, Li, and B-like g 3 ®Bi ions in 
fj,eW. The total nuclear-polarization contrihution to the specific differ¬ 
ences Ahl^'Hfs.np and AhbLIhfs.np are also presented in /leV. 


State 

Is 

2 s 

to 

AblF'HFS.NP Aggi^HFS.NP 

Ai^HFS.NP 

50.34 

55“ 

8.865 

0.730 


Ai?HFS,SNP 


-0.340 

-0.095 

0.025(12) -0.09(5) 


“Ref. (H. 


screened NP corrections Ai^HFS.NP and A£’hfs,Snp, respec¬ 
tively, are presented in Table |III] Here, we employ the nu¬ 
clear single-particle model for the description of the Bohr- 
Weisskopf effect. In the case of the Is state a fair agree¬ 
ment has been achieved with the previous value ifl^ . that was 
obtained within the point-like magnetic moment approxima¬ 
tion. Now we introduce two specific differences: the first 
is between the HFS of H- and Li-like ions AJjlE'hfs de¬ 
fined as AJjl^'HFS = ^'E'hfs^ “ ChlA i^j^pg, and the sec¬ 
ond is between the HFS of H- and B-like ions A^gFinFS = 

_ ^ggAFigpg. The parameters ^hl and 
^hb are chosen in a way to cancel the Bohr-Weisskopf ef¬ 
fects, and they appear to be rather stable with respect to vari¬ 
ations of the nuclear model of the magnetization distribution 
(0. In the case of | 3 ®Bi these parameters are chosen to be 
?HL = 0.16886 and <^hb = 0.014459. In Table Hill we present 
the obtained results for the total NP contribution to the spe¬ 
cific differences under consideration. As we can stress now 
the NP effects are essentially reduced in the both differences. 
The obtained results for the AgLFiHFS.NP and Aggi^HFS.NP 
appear to be rather stable with respect to the changes of the 
charge and magnetization distribution models. In view of this 
we assign an uncertainty of 50% to the total NP corrections 
to the specific differences. This result opens the possibility 
for further theoretical improvements of the HFS specific dif¬ 
ferences and tests of the magnetic sector of strong-held QED. 
Moreover, this can lead to the determination of the nuclear 
magnetic moments from a comparison between the theoreti¬ 
cal and experimental values for the specific differences. 


We can now also consider the situation for the bound- 
electron g factor. The leading order and the screened NP 
corrections are given by the same diagrams as for the HFS, 
see. Figs. [3 and |4] respectively. In Table |IV] we present 
our numerical results obtained for the g factors of H-, Li-, 
and B-like lead ions. Our values for AgNP are in a reason¬ 
able agreement with the previous calculations of Ref. MlSll . 
We also present the results for the specific differences be¬ 
tween the H- and Li-like g factors ggg defined by ggg = 
g[(is)^ 2 s] _ g[^®l, and between the H- and B-like g fac¬ 


tors ggg = g[(i®)^( 2 s)^ 2 pi/ 2 ] _ ^ggg[i®l. In the case of Pb 
ions the ^ parameters are chosen to be ^hl = 0.1670264 M 
and $hb = 0.0097416 fill] . As one can see from the table the 
NP corrections are canceled in a specific differences by about 


TABLE IV: Individual contributions to the leading order AgNP and 
the screened AgsNP nuclear-polarization corrections to the ground 
state g factor of H-, Li, and B-like 82 ®Pb ions in 10“®. The 
total nuclear-polarization contrihution to the specific differences 
AgHL.NP snd AgHB NP ^6 slso presented in 10~®. The leading 
order nuclear-polarization corrections are compared with the previ¬ 
ous calculations d- 


State 

Is 

2 s 

2 Pl/2 

AgHL.NP 

AgHB.NP 

AgNP 

-19.77 

-22“ 

-3.444 

-3.8“ 

-0.291 

-0.32“ 



AgSNP 


0.129 

0.104 

-0.013(6) 

0.006(3) 

0.004(6)*’ 


“Ref, ^Ref. fll. 


two orders of magnitude. Our conservative estimation of the 
total uncertainty is of about 50% of the effect. In compari¬ 
son with the rough estimation of the screened NP correction 
made in Ref. Il111 , here we rigorously take into account the 
first-order interelectronic-interaction correction to the nuclear 
polarization. We have also evaluated the effect originating 
from nuclear excitations induced via via two types of mag¬ 
netic interactions: one is due to the interaction with a constant 
magnetic field, and the other one is due to the magnetic in¬ 
teraction with an electron. In Ref. iH such corrections were 
referred to as nuclear magnetic susceptibility corrections to 
the g factor. The contributions of this effect to the g factor are 
found to be —6.6 x 10“^^, —1.1 x 10“^^, and —0.4 x 10“^^ 
for the Is, 2s, and 2pi/2 states, respectively. The correspond¬ 
ing contributions to the specific differences are at least by an 
order of magnitude smaller than our uncertainty. Thus, we 
can state an improvement of the accuracy of NP correction 
to the specific difference ggg by a factor of 2. In view of 
this new result we can push back the ultimate limit, defined 
by the NP effects, to the specific difference ggg and conse¬ 
quently to the possible accuracy of the determination of the 
fine structure constant. The recommended value of the fine 
structure constant according to the recent CODATA 13211 is 
a = 1/137.035 999 074(44). The corresponding uncertainty 
in the specific difference <5ggg[a] = 5 x 10“^^ is 1.5 times 
larger than the theoretical limit given by the NP uncertainty 
i5ggB NP = 3 X 10“^^. Another principal uncertainty in the 
specific difference is coming from the remaining finite nuclear 
size effect. However, this uncertainty can be substantially re¬ 
duced by employing the more accurate charge distribution pa¬ 
rameters obtained from muonic atoms , 34 ]. 


To conclude, we have evaluated the leading order and the 
screened NP corrections to the binding energies, HFS, and 
bound-electron g factor of heavy highly charged ions. The 
interelectronic-interaction effects have been rigorously evalu¬ 
ated within the QED perturbation theory up to the first order 
in 1/Z. The effect of the nuclear polarization has been evalu¬ 
ated for the specific differences constructed in a way to cancel 
the nuclear size corrections. In all cases considered here it 
turns out, that the NP corrections determining the ultimate ac¬ 
curacy cancel substantially. Therefore, the rigorous investiga- 
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tions of the specific differences provide a unique opportunity 
to test the strong-held QED with much higher accuracy than 
expected before. The ultimate accuracy of the nuclear polar¬ 
ization for the specihc difference has been now improved 

by a factor of 2. This may clear the way for more accurate 
determination of the hne structure constant from the strong- 
held QED with a precision similar to the one obtained from 
the investigations of the free-electron g factor. 
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